Generalised linear mixed model analysis via sequential Monte Carlo sampling
نویسنده
چکیده
We present a sequential Monte Carlo sampler algorithm for the Bayesian analysis of generalised linear mixed models (GLMMs). These models support a variety of interesting regression-type analyses, but performing inference is often extremely difficult, even when using the Bayesian approach combined with Markov chain Monte Carlo (MCMC). The Sequential Monte Carlo sampler (SMC) is a new and general method for producing samples from posterior distributions. In this article we demonstrate use of the SMC method for performing inference for GLMMs. We demonstrate the effectiveness of the method on both simulated and real data, and find that sequential Monte Carlo is a competitive alternative to the available MCMC techniques.
منابع مشابه
A sequential Monte Carlo approach for MLE in a plant growth model
Parametric identification of plant growth models formalized as discrete dynamical systems is a challenging problem due to specific data acquisition (system observation is generally done with destructive measurements), non-linear dynamics, model uncertainties and highdimensional parameter space. In this study, we present a novel idea of modeling plant growth in the framework of non-homogeneous h...
متن کاملParsimonious classification via generalised linear mixed models
We devise a classification algorithm based on generalised linear mixed model (GLMM) technology. The algorithm incorporates spline smoothing, additive model-type structures and model selection. For reasons of speed we employ the Laplace approximation, rather than Monte Carlo methods. Tests on real and simulated data show the algorithm to have good classification performance. Moreover, the result...
متن کاملBayesian Inference for Generalized Additive Regression based on Dynamic Models
We present a general approach for Bayesian inference via Markov chain Monte Carlo MCMC simulation in generalized additive semiparametric and mixed models It is particularly appropriate for discrete and other fundamentally non Gaussian responses where Gibbs sampling techniques developed for Gaussian models cannot be applied We use the close relation between nonparametric regression and dynamic o...
متن کاملDefault Bayesian model determination methods for generalised linear mixed models
A default strategy for fully Bayesian model determination for GLMMs is considered which addresses the two key issues of default prior specification and computation. In particular, the concept of unit information priors is extended to the parameters of a GLMM. A combination of MCMC and Laplace approximations is used to compute approximations to the posterior model probabilities to find a subset ...
متن کامل